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The statistics of extrasolar planetary systems indicate that the de¬ 
fault mode of planet formation generates planets with orbital peri¬ 
ods shorter than 100 days, and masses substantially exceeding that 
of the Earth. When viewed in this context, the Solar System is 
unusual. Here, we present simulations which show that a popular 
formation scenario for Jupiter and Saturn, in which Jupiter migrates 
inward from a > 5 AU to a ~ 1.5 AU before reversing direction, can 
explain the low overall mass of the Solar System's terrestrial plan¬ 
ets, as well as the absence of planets with a < 0.4 AU. Jupiter's 
inward migration entrained s > 10 — 100 km planetesimals into low- 
order mean-motion resonances, shepherding and exciting their orbits. 
The resulting collisional cascade generated a planetesimal disk that, 
evolving under gas drag, would have driven any pre-existing short- 
period planets into the Sun. In this scenario, the Solar System's ter¬ 
restrial planets formed from gas-starved mass-depleted debris that 
remained after the primary period of dynamical evolution. 

planetary dynamics | solar system | extrasolar planets 


Significance 

T he Solar System is an unusual member of the galactic 
planetary census in that it lacks planets that reside in 
close proximity to the Sun. In this work, we propose that the 
primordial nebula-driven process responsible for retention of 
Jupiter and Saturn at large orbital radii and sculpting Mars’ 
low mass is also responsible for clearing out the Solar Sys¬ 
tem’s innermost region. Cumulatively, our results place the 
Solar System and the mechanisms that shaped its unique or¬ 
bital architecture into a broader, extrasolar context. 


Introduction 

A full understanding of the formation and the early evolu¬ 
tion of the Solar System ranks among natural science’s grand 
challenges, and at present, even the dominant processes re¬ 
sponsible for generating the observed planetary architecture 
remain elusive. Nonetheless, the past three decades have gen¬ 
erated remarkable progress o, and critically, the discovery of 
thousands of extrasolar planets has placed the Earth and the 
Solar System into the broader context of the galactic planetary 
census. 

Perhaps the most important exoplanet-related discovery 
has been the realization that roughly half of the Sun-like 
stars in the solar neighborhood are accompanied by systems 
of one or more planets on low-eccentricity orbits with peri¬ 
ods ranging from days to months, and masses falling in the 
1 Me < Mp < 50 Me range ([2l [^. This dominant popula¬ 
tion of planets (which often presents tightly packed, nearly 
co-planar multiple systems) contrasts sharply with the Solar 
System, whose inner edge is marked by Mercury’s 88-day (0.4 
AU) orbit (see Figure 1). An iconic example from the new 
planetary catalog is the Kepler-11 system, which encompasses 
at least six planets comprising more than ~ 40 Earth masses 
In short, the exoplanetary surveys have revealed a hith¬ 
erto unrecognized oddity of the Solar System. Relative to 



Fig. 1. Orbital distribution of sub-Jovian extrasolar planets. A 
collection of transiting planet candidates with radii R < 5/^0, de¬ 
tected by the Kepler mission is shown. The radial distance away 
from the center of the figure represents a logarithmic measure of 
the planetary semi-major axis, such that the origin corresponds to 
the Sun’s surface. The sizes of the individual points represent the 
physical radii of the planets. Further, the points are color-coded 
in accordance with multiplicity. The orbits of the terrestrial plan¬ 
ets are also shown. Despite observational biases inherent to the 
observed distribution (e.g. transit probability, detectability) that 
work against detection of planets at increasing orbital radii, the 
raw contrast to our own Solar System is striking. 


other Sun-like, planet-bearing stars, our terrestrial region is 
severely depleted in mass. 

Reserved for Publication Footnotes 


^See refs miD for an alternative view. 

notable system within the resonant extrasolar population is GJ 876, where the inner planet is 
substantially less massive than the outer. In accordance with the picture of resonant transport de¬ 
lineated in ref. (T8\ . this system likely failed to satisfy the conditions required for migration reversal 
and decayed to a compact orbital configuration fT^ . 
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A few related peculiarities are also evident within the inner 
Solar System. Specifically, cosmo-chemical evidence suggests 
that while the fundamental planetary building blocks (plan- 
etesimals) formed within ~ 1 Myr of the Sun’s birth ([5}, the 
final assembly of the terrestrial planets occurred on a timescale 
of 100 — 200 Myr, well after the dispersal of the nebular gas 
This is at odds with the inferred compositions of extrasolar 
Super-Earths, which are thought to have substantial gaseous 
atmospheres. Additionally, the exceptionally small masses of 
Mercury and Mars sugges^that the terrestrial planets formed 
out of a narrow annulus of rocky debris, spanning 0.7—1 AU 
(where lAU is the mean distance between the Earth and the 
Sun) ([9} . Such a narrow annulus is at odds with so-called 
minimum mass Solar Nebula ([TOlfTT]) . 

Within the framework of a radially confined solid component 
of the inner Solar nebula, the inner edge of the annulus is en¬ 
tirely artificial. Indeed, at present there exists no compelling 
justification for its origin. A plausible explanation may stem 
from the dynamical evacuation of solid material by a popula¬ 
tion of primordial close-in planets cai- We shall investigate 
this possibility further in this study. 

Unlike the inner edge of the annulus, a body of previous 
work has demonstrated that the outer edge can be naturally 
sculpted by inward-then-outward migration of Jupiter m- 
Within protoplanetary disks, long-range migration of giant 
planets results from tidal interactions with the nebula and 
viscous transport (na. Eor single planets, orbital evolution 
is typically inwards. However, the process of resonant lock¬ 
ing between two convergently migrating planets can lead to a 
reversal of the migration direction dm. 

The process of resonant migration reversal for gap-opening 
planets (i.e. objects with M > Mjup) is a well-understood 
result of planet-disk interactions, and only requires the outer 
planet to be somewhat less massive than the inner. To this 
end, it is worth noting that all of the known mean-motion com¬ 
mensurate pairs of giant planets that reside beyond a > 1 AU 
have the more massive object on the inside dnilTel, suggesting 
that the operation of this mechanism is widespreacQ 

Within the Solar System, it is inferred that Jupiter initially 
migrated inwards from its primordial formation site (presum¬ 
ably 3 — 10 AU) to ~ 1.5 AU, and subsequently reversed its 
evolutionary track as a consequence of locking into a 3:2 mean- 
motion resonance with a newly-formed Saturn. This special 
case of the generic resonant migration reversal mechanism is 
informally referred to as the “Grand Tack” scenario (|13|) . In 
addition to the aforementioned truncation of the inner solid 
nebula, this putative sequence of events is attractive in that 
it naturally explains how the Solar System’s giant planets 
avoided spiraling into Sun m, accounts for the origins of 
compositional differences within the Asteroid belt (O, pro¬ 
vides a mechanism for delivery of water into the terrestrial 
region (HSl, and generates a compact orbital configuration 
needed for the subsequent instability-driven orbital evolution 
of the outer Solar System dm HD). 


Resonant Transport and Collisional Evolution 

An early inward migration for Jupiter has a number of reper¬ 
cussions that come to light when one places the Solar System 
into the broader context provided by the observations of extra¬ 
solar planets. An inescapable consequence of Jupiter’s trek is 
the resonant capture and the ensuing inward entrainment and 
transport of solid material (|23D . In particular, when a given 
planetesimal’s orbital period becomes nearly rational multiple 
of the orbital period of Jupiter, gravitational perturbations 
become coherent and force the planetesimal to maintain the 
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Fig. 2. Orbital evolution of planetesimals embedded in the So¬ 
lar nebula, under the effects of a migrating Jupiter. As Jupiter 
moves inwards from 6 AU to 1.5 AU, planetesimals are swept up by 
mean motion resonances (MMRs). Panel A of this figure shows the 
increase in the planetesimal eccentricity associated with resonant 
transport. Note that at the end of the Jupiter’s trek, there exists 
a strong enhancement in the planetesimal density at the Jovian 2:1 
MMR. Panel B depicts the preferential population of Jupiter’s in¬ 
terior MMRs. Each planetesimal in the simulation is color-coded 
in accord with its initial condition, and the resultant curves track 
the orbital excursions of the small bodies as Jupiter’s orbit shrinks. 
Jupiter’s return to ~ 5 AU is not modeled directly. In the presented 
simulation, we assumed a planetesimal size of s = 100 km. Similar 
figures corresponding to s = 10 km and s = 1000 km can be found 
in the SI. 


same period ratio, leading to a decrease in the planetesimal 
orbit’s semi-major axis dm. The most common commensura- 
bility at which capture occurs is 2:1, although numerous other 
possibilities exist. 

In order for resonant interactions to be effective, the plan¬ 
etesimal in question must not be hydrodynamically coupled 
to the nebular gas on the orbital timescale. In practice, this 
means that the planetesimal size must significantly exceed 
s ^ 1 cm ([25l) . Modern calculations of planetesimal formation 
suggest that planetesimals grow very rapidly to radii of order 
s ~ 100 km (|26D and do not experience severe interactions 
with the nebula m- Under such conditions, the entire solid 
component of the nebula swept up by Jupiter’s resonances 
will be captured and transported inwards. In particular, if 
one adopts the conservative assumption of a minimum-mass 
Solar nebula, and an oft-quoted solid-to-gas ratio of ~ 1%, the 


^See SI for similar calculations corresponding to s = 10 km and 1000 km planetesimals. 
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total amount of mass swept up by the resonances is ~ 10 — 20 
Earth masses, about an order of magnitude greater than the 
cumulative mass of the terrestrial planets. 

A planetesimal that is driven inward by resonant migra¬ 
tion experiences a concomitant increase in its orbital eccentric¬ 
ity. This effect is generally well understood and stems from 
adiabatic invariance (see Supporting Information (SI) for a 
discussion). The eccentricity of a planetesimal embedded in 
a gaseous nebula cannot grow indefinitely however, as it is 
damped by aerodynamic drag. Thus, the adiabatic excita¬ 
tion of eccentricity typically stalls at a finite value dependent 
on the planetesimal size (which controls the magnitude of the 
dissipative forcing). 

We have calculated the orbital evolution of s = 10 km, 100 
km, and 1000 km planetesimals as they are swept up by an 
inward-migrating Jupiter. The results of calculations corre¬ 
sponding to the nominal case of s = 100 km are presentecQin 
Figure (2). The simulations suggest that for the aforemen¬ 
tioned planetesimal sizes, the equilibrium eccentricities are 
e ~ 0.2, 0.3 and 0.5 respectively. Therefore, substantial or¬ 
bital crossing will occur between planetesimals caught in reso¬ 
nances with Jupiter and those residing within as-yet unswept 
regions of the disk. 

The simulations reported in this work were terminated upon 
Jupiter’s arrival to 1.5 AU. As already mentioned earlier, a res¬ 
onant encounter with Saturn followed by reversal of migration 
is envisaged to have occurred subsequently. However, outward 
migration is not important to the problem at hand because in¬ 
terior material cannot be transported to wider orbits by this 
process. 

Solid bodies on crossing orbits within densely populated 
disks experience collisions, which can result in either accre¬ 
tion or fragmentation. The outcome is principally determined 
by the specific energy of the impact: if this quantity exceeds 
a critical value characteristic of catastrophic disruption, the 
target is shattered into two or more pieces ([^ [2^ . Adopt¬ 
ing parameters appropriate for high-velocity impacts among 
strong basaltic objects in the gravity-dominated regime (see 
SI for these parameters), we find analytically that the spe¬ 
cific impact energy safely exceeds its threshold value across 
the range of planetesimal sizes invoked above for impactor-to- 
target mass ratio of ~ 0.1 or greater (see SI). In other words, 
our results suggest that even though one may expect that the 
real planetesimal disk will harbor a distribution of planetes¬ 
imal sizes, the trade-off between size-dependent orbital exci¬ 
tation and threshold impact energy leads to an environment 
where objects of any size above s > 10 km can be destroyed 
by bodies that are ~ 10 times less massive. Jupiter’s reso¬ 
nant shepherding of planetesimals thus initiates a collisional 
cascade (TO that grinds down the planetesimal population to 
smaller sizes. 

Although the details of resonantly-forced collisional grind¬ 
ing can be complex, an important feature of this process is 
that once the size of a given planetesimal population is dimin¬ 
ished to a point where the effects of aerodynamic drag become 
important (for example, s < 1 km at 1 AU), the planetesimals 
will experience a runaway inwards drift m- Importantly, the 
same process facilitates the removal of material from Jovian 
resonances and thereby yields a critical planetesimal size be¬ 
low which collisional grinding subsides. Thus, the collisional 
cascade initiates well before Jupiter reaches its innermost tack, 
and proceeds as long as Jupiter’s migration direction is inward. 

Given the exceptionally large impact frequency expected 
within a mature protoplanetary disk and the dominantly de¬ 
structive nature of collisions discussed above, we expect that 
a sizable fraction, if not all, of the transported population of 
planetesimals will be disrupted and undergo rapid orbital de¬ 
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Fig. 3. Orbital decay of a hypothetical compact system of Super- 
Earths (an analog of the Kepler-11 system) residing within the ter¬ 
restrial region of the primordial Solar System. Following a colli¬ 
sional avalanche facilitated by Jupiter’s migration, a population of 
planetesimals (here assumed to be ground down to s = 100 m) de¬ 
cays inwards and resonantly shepherds the interior planets into the 
star. Planetesimal orbits are shown with colored lines, while the 
planetary orbits are shown with black and gray lines. Specifically, 
the planetary semi-major axes are shown in black while the perihe¬ 
lion and aphelion distances are shown in gray. Note that the results 
shown herein are largely independent of planetesimal size, as long as 
the planetesimals are small enough to drift inwards on a timescale 
smaller than ~ 1 Myr due to aerodynamic drag. 

cay following Jupiter’s reversal of migration direction. This 
feature is of critical importance for explaining the Solar Sys¬ 
tem’s lack of close-in Super-Earths. 


Decay of Primordial Close-in Planets 

The dominant formation channel (distant formation followed 
by extensive inward migration vs. in-situ conglomeration 
l|32l |33l) for extrasolar Super-Earths remains controversial. 
However, a generally agreed-upon framework of core-nucleated 
accretion of giant planets dictates that the formation of solid 
multi-Earth-mass cores precedes the formation of giant plan¬ 
ets da. Thus, given that the formation of tightly packed 
close-in systems is ubiquitous in the galaxy, it can be reason¬ 
ably speculated that at the time of Jupiter’s inward journey, 
a similar population of first-generation planets existed in the 
Solar System. If planets formed from such material, however, 
they were destroyed. 

In exactly the same way as an inward-migrating Jupiter 
captures planetesimals into resonance, inward-migrating plan¬ 
etesimals will lock into resonance with close-in planets. Pro¬ 
vided that the cumulative mass of the resonant planetesimal 
population is not negligible compared to the mass of the close- 
in planets, the planetesimals will gravitationally shepherd the 
close-in planets into the Sun. In other words, the inward-then- 
outward migration of Jupiter in the early Solar System wiped 
the inner Solar System’s slate clean, setting the stage for the 
formation of a mass-depleted, gas-starved second generation 
of terrestrial planets (Hnsi. Indeed, within the framework of 
this picture, the material from which the Earth formed is ei¬ 
ther the remainder of the violent collisional avalanche, or has 
been largely emplaced by Jupiter’s outward migration. 


^This example is used for definitiveness. We are not suggesting that a primordial population of 
the Solar System’s close-in planets would have necessarily borne any similarity to the Kepler —11 
system. 
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To illustrate the above process, we examined the dynamical 
evolution of the Kepler—11 planetary systenj^ when placed 
within the inner edge of the Solar nebula, and under the 
gravitational influence of an extensive population of exterior, 
inward-drifting planetesimals. The computed evolutionary se¬ 
quence is shown in Figure (3). Clearly, dissipative resonant 
transport provides an efficient mechanism for driving close-in 
planets into the central star. Indeed, the sequence of events as¬ 
sociated with Jupiter’s so-called “Grand Tack” may well have 
constituted a veritable grand attack on the Solar System’s 
original population of short-period Super-Earths. 


Discussion 

This scenario provides a natural explanation for why the inner 
Solar System bears scant resemblance to the ubiquitous multi¬ 
planet systems discovered by the Doppler velocity surveys and 
by the Kepler mission. Moreover, the physical processes that 
we invoke (namely giant planet migration, collisional disrup¬ 
tion of planetesimals, aerodynamic drag, and resonant shep¬ 
herding) are generic. In consequence, the mechanism de¬ 
scribed herein should also operate within a non-negligible frac¬ 
tion of extrasolar planetary systems. Accordingly, a series of 
observational predictions can be formulated. 

First, our calculations imply a strong anti-correlation be¬ 
tween the existence of multiple close-in planets and giant plan¬ 
ets at orbital periods exceeding ~ 100 days within the same 
system. The existing exoplanet catalog is not yet sufficiently 
detailed to test this hypothesis (HU. However, direct assess¬ 
ment of the validity of this prediction will be provided by the 
upcoming TESS and K2 missions. Second, the spectral en¬ 
ergy distributions of protoplanetary disks hosting gap-opening 
planets should exhibit strong infra-red enhancements (ISH , as 
a consequence of collisional heating and the associated pro¬ 
duction of dust. Moreover, dust emission morphologies in 
such disks could in principle exhibit asymmetrical structure 
([36|). Most dramatically, our work implies that the majority of 
Earth-mass planets are strongly enriched in volatile elements 
and are uninhabitable. 
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Supplementary Information 


Methodology 

Numerical simulations reported in this work made use of a 
modified mercuryO N-body integration software package ([T]). 
The standard gravitational dynamics solver was augmented 
with the auxilirary effects of aerodynamic drag (HEli for small 
objects as well as fictitious forces that mimic disk-driven mi¬ 
gration and circularization of planets. The hybrid symplectic 
Bulirsch-Stoer algorithm @ was utilized throughout. 


Inward Migration of Jupiter 

Because of dynamical evolution that occurred during the 
nebular epoch of the Solar System’s lifetime, Jupiter’s current 
orbital radius is not informative of its formation site. How¬ 
ever, theoretical arguments suggest that the formation of giant 
planets generally occurs beyond the snow-line i.e. at a distance 
of ~ a few astronomical units dsiisi. Accordingly, in this work 
we take Jupiter’s initial semi-major axis to be (aj)i = 6 AU. 

We employ the standard assumption that the viscosity of 
the Solar nebula was not overwhelmingly large, meaning that 
upon formation Jupiter opened a substantial gap in the pro¬ 
toplanetary disk and subsequently migrated inwards in the 
so-called “type-II” regime ©• Because type-II migration is 
associated with the global evolution of the disk, the charac¬ 
teristic migration timescale 


2 1 da 

Ta a dt 


[ 1 ] 


is synonymous with the viscous timescale of the disk. Thus, 
in agreement with published hydrodynamical simulations, we 
adopted a value of (Ta)j = 2 x 10^ years ([8}, corresponding to 
a Shakrua-Sunayev a ~ 0.003. 

To facilitate migration, we introduced an acceleration amig 
of the form 


a-mig — , [2] 

Ta 

where v is the orbital velocity. The simulations were termi¬ 
nated when Jupiter reached its envisioned tacking semi-major 
axis of (aj)f = 1.5 AU. 

Planetesimal Evolution 

As shown in Figure (2) of the main text, a dominant frac¬ 
tion of the planetesimal swarm swept by Jupiter’s interior 
mean-motion commensurabilities was captured into the cor¬ 
responding resonances. The conditions for resonant capture 
of small objects initially residing interior to Jupiter’s orbit 
can be deduced analytically within the framework of the cir¬ 
cular restricted three body problem (nnj. Ignoring the effects 
of Jupiter’s minuscule eccentricity, adiabatic theory dictates 
that capture into a k : k — 1 resonance is guaranteed if the 
initial planetesimal eccentricity is less than a critical value: 

!=>! 

where (^k is a interaction coefficient that although different for 
each resonance, is typically of order unity (for interior 2:1 and 
3:2 resonances, ^2 = —1.19 and ^3 = —2.02 respectively) (pTl) . 
It can be easily checked that provided reasonable parameter 
choices, this conditior0is almost certainly satisfied in the pri¬ 
mordial Solar nebula. 

An additional requirement for capture to be certain is the 
fulfillment of the adiabatic condition. For the problem at 
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Fig. 4. Evolutionary sequence of s = 10 km planetesimals. This 
figure is an analog of Figure (2) of the main text. 


hand, this means that the orbital convergence timescale must 
be longer than the resonant libration period, defined as ca: 


Piib = — 
n 


k-1 


4/3 


Mj 

Mo 


Cl 




TT- CfcC 


1 - 1/2 


[4] 


where e explicitly refers to the eccentricity at the libration 
center (i.e. resonant equilibrium). 

The first resonant trajectories (associated with the emer¬ 
gence of a homoclinic curve in phase space) appear at an equi¬ 
librium eccentricity of 


cl 


1/9 


\Mq) ■ 


The corresponding resonance width is given by 


Aa 


8 / (fc-i)ci y^' (^\ 

32/3 k^ ) \Mq) 


2/3 


[5] 

[ 6 ] 


Accordingly, setting the dissipative resonance crossing time 
equal to the libration period, we obtain the criterion for adia- 


^Quantitatively, the critical eccentricity evaluates to e^rit — 0.15 and 0.12 for 2:1 and 3:2 
resonances respectively. 

®Here, we assume a starting eccentricity of e = 0. 
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Fig. 5. Evolutionary sequence of s = 1000 km planetesimals. 
This figure is an analog of Figure (2) of the main text. 


batic evolution: 

^ > _ 27r/nj _ / (kMj \ 

~ 4(3)2/3(yfc(fc - l))2/9 \ Me J 




Quantitatively, the above expression dictates that in order for 
resonant capture to be ubiquitously successful, the migration 
timescale must exceed Ta > 10^ years. Indeed, the typically 
quoted range of type-II migration timescales is in agreement 
with the adiabatic condition cni. 

Small bodies trapped in resonances will maintain a con¬ 
stant period ratio and will therefore migrate inwards along 
with Jupiter. Associated with this migration is an adiabatic 
invariant dlsllll 


^/a — {k — l)\/l — e^j = const. [ 8 ] 

The conservation of this quantity implies that as the particle 
semi-major axis is decreased, its eccentricity must grow. 

By taking a derivative of the above expression, we can obtain 
a differential equation for test particle’s eccentricity growth. 
Combined with the requirement that the planetesimal semi¬ 
major axis must decrease in concert with that of Jupiter [see 
equation once the resonance is established, the differential 
equation can be solvecQto obtain the planetesimal eccentricity 
as a function of time: 


This expression implies that the characteristic timescale for 
any given particle to reach its equilibrium eccentricity is of 
order ~ 0.1 — 0.5Ta. 

Within the context of our numerical simulations, the solid 
component of the nebula was initialized as a disk comprised of 
1000 objects with inner and outer radii of (adisk)in = 0.1 AU 
and (adisk)out = 6 AU respectively. The outer radius of this 
disk is not indicative of the true outer radius of the Solar neb¬ 
ula. Rather, planetesimals with semi-major axes greater than 
Jupiter’s initial semi-major axis were ignored in our simula¬ 
tions because they cannot be transported inwards by sweeping 
exterior mean motion resonances (divergent resonant encoun¬ 
ters cannot lead to capture) fTOl ITT]) . 

The radial distribution of the planetesimals followed a sur¬ 
face density profile equivalent to that of a Mestel disk (O: 

S = 110 ] 

where Uq is the disk surface density at a semi-major axis ao. 
We note that the total planetesimal mass transported inwards 
by resonant shepherding (assuming that first-order resonances 
such as 2:1 and 3:2 are primarily responsible for long-range mi¬ 
gration) is 


Mtot — 


J(a.), ^ CL J 


a da 


' («j)f 

27r e / Eo ao (aj)72^'’^, 


[ 11 ] 


where / ~ 0.01 is the disk solid-to-gas mass ratio and e < 1 is 
the resonant transport efficienc^Q Quantitatively, only ~ 20% 
of the captured objects were observed to chaotically escape the 
resonances during the entirety of Jupiter’s inward trek in the 
simulations. 

Adopting MMSN-like parameters of So — 2000 g/cm^ at 
ao = lAU, we obtain Mtot ~ 2 OM 0 ^ Mj. Consequently, 
the gravitational back-reaction of the planetesimals onto an 
inward-migrating Jupiter is almost certainly unimportant, and 
for the purposes of this set of simulations we ignored this ef¬ 
fect. Correspondingly, Jupiter’s eccentricity remained at a 
near-null value throughout the integrations. 

Assuming that the solid and gaseous component of the neb¬ 
ula follow the same surface density profile , radial pressure 
support will yield a sub-Keplerian circular velocity profile of 
the form (E}: 


Vgas = VK\ 1 - 3^ (f = Vk{ 1 - rf) ip, 


[ 12 ] 


where Cg and ^'K are the sound and Keplerian speeds respec¬ 
tively, while p) denotes the azimuthal unit vector. Prototypical 
disk parameters yield a value of 77 in the range 0.001 — 0 . 01 . 
Accordingly, in this work we adopted 77 = 0.005 (correspond¬ 
ing to a disk aspect ratio of h/r ~ 0.05) for all calculations. 

Planetesimals embedded in a sub-Keplerian gas disk will 
experience orbital decay due to aerodynamic drag. The hy¬ 
drodynamic Reynolds number relevant for planetesimals big¬ 
ger than s > 10 m substantially exceeds unity dm. In this 
regime, drag acceleration adrag is quadratic in velocity and 
independent of fluid viscosity Cl: 


ttCd 2 

^drag — T S Pgas'^rel Vrel • 

2m 


[13] 


1 — exp(i/ra) — 2fc(l - exp(i/(2Ta)) 


(fc-l)2 


[9] 


^Note that e is primarily limited by chaotic diffusion and the associated removal of objects from 
Jovian resonances fT^ . 
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In the above expression, Vrei = v — Vgas is the relative veloc¬ 
ity between the planetesimal and the gas, whereas m is the 
planetesimal mass. For the entire set of our simulations, the 
aforementioned acceleration was implemented with a drag co¬ 
efficient of Cd — 0.5. 

In addition to the decrease in semi-major axes, the introduc¬ 
tion of this acceleration results in circularization of the orbits 
and damping of the vertical motion (inclinations). To lead¬ 
ing order in e and z, the rates of these dissipative effects are 
expressed as follows 


where 


1 da 
a dt 
1 de 
e dt 


2 /5 2^ 1-2^ 2 

= , 1^2 
i dt To V 8 2 




To = 


TtCd 2 

S Pga 


2m 


,vk 


[14] 


[15] 


In light of equations [^, it is important to recall that al¬ 
though first order resonant interactions with Jupiter will mod¬ 
ulate the planetesimal eccentricity [equation 13, they will not 
excite orbital inclination (HD- Aerodynamic drag on the other 
hand, tends to damp vertical motion of eccentric orbits. Thus, 
the process of resonant transport leads to a dissipative confine¬ 
ment of the solid component of the nebula to the mid-plane 
of the disk. 

Note that here, we are assuming that turbulent fluctuations 
in the disk are not sufficiently large to obstruct resonant cap¬ 
ture. This is likely a safe assumption given that the orbital 
plane of planetesimals will likely be confined to the dominantly 
laminar mid-plane. However, an improved iteration of the con¬ 
sidered model could in principle account for this effect directly. 

A promising theory for the dominant channel for planetesi¬ 
mal formation invokes the gravitational collapse of solid over¬ 
densities within the disk GSl, which are in turn fabricated by 
turbulent forcing and streaming instabilities ([20|. Within the 
framework of this picture, planetesimals are born big, with 
peak masses of order that of Ceres. As a typical planetesimal 
size, it is sensible to choose s = 100 km and simulations em¬ 
ploying this default value are reported in the main text. For 
completeness, we have repeated the numerical experiments of 
inward migration of Jupiter and the associated transport of 
planetesimals with s = 10 km and s = 1000 km. The corre¬ 
sponding results are shown in Figures (4) and (5) respectively. 

As expected, hydrodynamical damping of smaller objects 
is substantially more pronounced, which means that resonant 
transport is more efficient at eccentricity excitation for higher- 
mass bodies. Accordingly, our simulations show that the equi¬ 
librium eccentricity attained by 10 km and 1000 km objects 
throughout their inward trek is e ~ 0.2 and e ~ 0.5 respec¬ 
tively. 


Collisions 

Resonant shepherding of planetesimals by Jupiter will lead 
to a concentration of small objects at particular (mean-motion 
resonant) orbits. In turn, such confinement will lead to a 
strong enhancement in the rate of collisions among planetesi¬ 
mals. A rough estimate of the collision frequency, v (defined 
as the inverse of the mean time between collisions experienced 
by a single particle), can be obtained from a simple nav cal¬ 
culation. Specifically, for the problem at hand we have (ED: 


— 


Mtot/m 

27r(e) tan(z)a^ 


7rs^z;K(e). 


[16] 
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Fig. 6. Orbital decay of a clone of the Kepler-11 system embed¬ 
ded within the inner Solar nebula. The calculation presented herein 
mirrors that shown in Figure (3) of the main text. However, within 
the framework of this simulation, dissipative interactions between 
the planets (but not planetesimals) and the disk are entirely ne¬ 
glected. Evidently, the effect of eccentricity damping on the planets 
plays a negligible role in determining the qualitative behavior of the 
system. 


If we take the entire planetesimal population to be com¬ 
prised of s = 100 km planetesimals, adopt (z) ~ 10“^ as sug¬ 
gested by the simulation^ and Mtot ~ IOM 0 , at a = I AU 
we obtain i/ ~ 0.05. This means that each member of the 
planetesimal population can be expected to suffer a collision 
each 20 orbits or so. This constitutes an exceptionally large 
collision rate and allows one to infer that collisional grinding 
of the small bodies initiated by Jupiter will be efficient. 

The detailed description of the outcome of collisions can in 
general be complex and may depend on numerous parameters 
that are specific to a particular impact. In an averaged sense, 
however, the controlling factor in determining the result of 
collisions is the specific energy ED: 



In this equation, m' is the impactor mass, M ^ m' is the 
target mass, and Venc is the encounter velocity. 

There exists a critical value characteristic of catas¬ 
trophic disruption, such that if Q exceeds QJ), the target is 
shattered into two or more pieces of mass no greater than M/2. 
In the gravity-dominated regime (appropriate for bodies with 
radius s > 0.5 km), the expression for reads ([231^ : 

Adopting p = 3 g/cm^, g = 0.5 and b — 1.36 appropriate for 
high-velocity impacts among basaltic objects and z;enc ~ ez;kep 
for the encounter velocity, we find that at lAU Q > for 
mass ratios of (m'/M) > 0.01,0.03 and O.I for s = 10,100 
and 1000 km respectively. Given that the process of colli¬ 
sional fragmentation acts to produce more objects of dimin¬ 
ished sizes, the above estimates suggest that the cascade en¬ 
abled by Jupiter’s inward migration will efficiently grind down 
the planetesimal population to sizes where orbital decay due 
to aerodynamic drag becomes an efficient removal mechanism. 


®Note that the smallest physically sensible value of the inclination that we can adopt for this 
calculation is fmin = tan~^ (s/a) ~ 10~®. 
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It should be noted that the simple calculations shown above 
ignore a number of detailed effects that can be taken into ac¬ 
count within the framework of a more sophisticated model. 
For example, we have ignored the various geometrical par¬ 
ticularities of collisions within resonant orbits, the associated 
collisional damping and viscous stirring, as well as the time- 
evolution of the planetesimal size distribution and the self- 
consistent removal of collisional fragments from Jovian reso¬ 
nances by rapid inward drift. While these issues may be of 
substantial interest as an avenue for follow-up work, they are 
not central to the arguments presented here, and may there¬ 
fore be omitted for the purposes of this study. 


Orbital Decay of a Kepler-11 Analog 

Extended inward flow of collisionally disrupted planetesi- 
mals may be halted (by resonant capture) if the migration 
route is obstructed by at least one massive planet m- How¬ 
ever, if the cumulative mass of planetesimals trapped in ex¬ 
terior planetary resonances becomes substantial (i.e. non- 
negligible with respect to the planetary mass), inward mi¬ 
gration of the entire system will occur. In light of this no¬ 
tion, we performed numerical experiments aimed at evaluat¬ 
ing the long-term fate of a compact, close-in aggregate of sub- 
Jovian planets embedded into the terrestrial region of the So¬ 
lar System, under the influence of an exterior flood of inward- 
spiraling planetesimals. 

Within the context of these simulations, the dynamical evo¬ 
lution of small bodies followed the same framework as that 
discussed above. The (presumed collisionally evolved) swarm 
of planetesimals was initialized as a circular annulus of 1000 
objects spanning 0.75 — 1.25 AU, and allowed to drift inwards. 
The total mass of the swarm was taken to be Mtot = 2 OM 0 . 

For the purposes of computing the effects of aerodynamic 
drag, a characteristic radius of s = 0.1 km was assumed for the 
entire population. We note that this value sets the timescale 
on which orbital decay takes place. In other words, similar 
dynamical evolution can be obtained by choosing a radius of 
s = 1 km and running the integrations for 10 times longer. 

As before, the maximal (aerodynamically-forced) orbital de¬ 
cay rates that satisfy the adiabatic condition can be calculated 
using classical perturbation theory CD- Assuming that the 
planetary orbit is circular, the libration period of a /c : /c — 1 
exterior resonance is m 


Piib = — 


n 



[19] 


where rh is the planetary mass, and (akin to Cfc) f,k is a co¬ 
efficient of order unity (for exterior 2:1 and 3:2 resonances, 
^2 = 0.428 and ^3 = 0.515 respectively). 

The equilibrium eccentricity at the inception of the reso¬ 
nance is 


^ ( ^kfh 

® 31/3 yk'^MQj 

with a corresponding resonance width of 


Aa 

- 8 1 

f ^kfh'\ 

a 

32/3/JI/3 1 

KMeJ 


2/3 


[ 20 ] 


[ 21 ] 


Calculating the resonance crossing time using equation p3] , 
we obtain the following expression for the adiabatic threshold: 


3Cd f (VKrf^ ^ - 1) ’^ f 


- 


- ^2/3 


, 4/3 


MqJ 


[221 


As an example, equation 22 dictates that adiabatic capture 


into a first order resonance with a m = IOM 0 planet residing 
in a MMSN at a = 0.4 AU is assured for particles with radii 
in excess of s > 1 km. While this is a generous constraint 
already, we note further that even in a trans-adiabatic regime, 
resonant capture will still occur, but with a diminished prob¬ 
ability m- To this end, our numerical simulations intention¬ 
ally break the adiabatic limit and thereby demonstrate that 
inward migration of close-in planets can still be successfully 
forced by debris substantially smaller than s < 1 km. 

In direct analogy with the capture of planetesimals into inte¬ 
rior Jovian resonances [equation]^, capture into exterior plan¬ 
etary resonances is assured if the particle eccentricity does not 
exceed a critical value (HU 


[23] 

\^k rn J 

Given that the hydrodynamic eccentricity damping timescale 
is ~ 100 times shorter than the semi-major axis decay 

timescale [see equations|14|, and the latter process is necessary 
for resonant capture in me first place, it is sensible to assume 
that this condition is well satisfied in the inner Solar nebula. 

As an emblematic example of a tightly packed set of close-in 
planets, we adopted the Kepler-11 system (j24p and inserted 
it interior to the decaying disk of planetesimals. The planets 
were initialized on their current near-circular orbits (specifi¬ 
cally, we adopted the “all-eccentric” orbital fit delineated in 
ref. ([24}). Gravitational interactions between the planets and 
the planetesimals were computed in the conventional N-body 
fashion. However, the self-gravity of the planetesimal disk was 
ignored to save computational costs. 

Unlike Jupiter, no fictitious migration due to the disk was 
imposed on the planets. Damping of the orbital eccentrici¬ 
ties due to angular momentum exchange with the disk was 
implemented via an acceleration of the form 


adamp = -2 *-^ [24] 

r^Te 

where as in equation [^, 

1 de 
e dt 

The quantity Te depends on numerous physical properties 
of the disk and is in general somewhat uncertain (HU). In this 
work, we performed simulations that adopted Te = 10 ^ years 
(shown in Figure 3 of the main text) as well as a run with 
Te = 00 (shown as Figure 6 ). Examination of the presented 
orbital solutions show that in either case, only mild eccen¬ 
tricities are attained as the planets are ushered towards the 
Sun. 

The exact extent of orbital excitation is set by the delicate 
interplay between adiabatic and dissipative dynamic^ How¬ 
ever, our results broadly suggest that for the problem at hand, 
the exact value of Te does not factor into the final answer ap¬ 
preciably. That is, our simulations show that irrespective of 
the details of the chaotic dynamical evolution that occurs as 
planetesimals lock into exterior resonance with the planets, 
the planetary orbits decay onto the host star over a timescale 
that is short compared with the typical lifetimes of protoplan¬ 
etary disks. 


®We note that for the non-restricted (i.e. planetary) three-body problem, an adiabatic invariant 
analogous to equation M can be defined fT^ . 


Footline Author 


PNAS I Issue Date | Volume | Issue Number | 9 






















’Grand'Attack” — 2015/4/2 — 0:19 — page 10 — #10 


References 

1. Chambers, J. E., A hybrid symplectic integrator that per¬ 
mits close encounters between massive bodies. Mon. Not. 
R. Astron. Soc. 304, 793-799 (1999) 

2. Weidenschilling, S. J. Aerodynamics of solid bodies in the 
solar nebula. Mon. Not. R. Astron. Soc. 180, 57-70 (1977) 

3. Adachi, L, Hayashi, C., Nakazawa, K. The gas drag ef¬ 
fect on the elliptical motion of a solid body in the pri¬ 
mordial solar nebula. Progress of Theoretical Physics 56, 
1756-1771 (1976) 

4. Press W. H. et ah. Numerical recipes in FORTRAN. The 
art of scientific computing. Cambridge: University Press, 
2nd ed. (1992) 

5. Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, 
J. J., Podolak, M., Greenzweig, Y. Formation of the Gi¬ 
ant Planets by Concurrent Accretion of Solids and Gas. 
Icarus 124, 62-85 (1996) 

6. Stevenson, D. J., Formation of the giant planets. Planet. 
Spac. Sci. 30, 755-764 (1982) 

7. Crida, A., Morbidelli, A., Cavity opening by a giant 
planet in a protoplanetary disc and effects on planetary 
migration. Mon. Not. R. Astron. Soc. 377, 1324-1336 
(2007) 

8. Walsh, K. J., Morbidelli, A., Raymond, S. N., O’Brien, 
D. P., Mandell, A. M. A low mass for Mars from Jupiter’s 
early gas-driven migration. Nature 475, 206-209 (2011) 

9. Papaloizou, J. C. B., Larwood, J. D., On the orbital evo¬ 
lution and growth of protoplanets embedded in a gaseous 
disc. Mon. Not. R. Astron. Soc. 315, 823-833 (2000) 

10. Henrard, J., apture into resonance - an extension of the 
use of adiabatic invariants. Cel. Mech. 27, 3-22 (1982) 

11. Murray, C. D., Dermott, S. F., Solar System Dynamics. 
UK: Cambridge University Press, (1999) 

12. Batygin, K., Morbidelli, A., Analytical treatment of plan¬ 
etary resonances. Astron. and Astrophys. 556, A28 (2013) 


13. Yu, Q., Tremaine, S. Resonant Capture by Inward- 
migrating Planets. Astron. J. 121, 1736-1740 (2001) 

14. Peale, S. J., Orbital resonances, unusual configurations 
and exotic rotation states among planetary satellites. 
Satellites 159-223 Tucson: University of Arizona Press 
(1986) 

15. Mestel, L., On the galactic law of rotation. Mon. Not. R. 
Astron. Soc. 126, 553 (1963) 

16. Wisdom, J., The resonance overlap criterion and the onset 
of stochastic behavior in the restricted three-body prob¬ 
lem. Astron. J. 85, 1122-1133 (1980) 

17. Malhotra, R. Orbital resonances in the solar nebula - 
Strengths and weaknesses. Icarus 106, 264 (1993) 

18. Landau, L. D., Lifshitz, E. M., Fluid mechanics. Course 
of theoretical physics Oxford: Pergamon Press (1959) 

20 , 

19. Johansen, A., Oishi, J. S., Mac Low, M.-M., Klahr, H., 
Henning, T., Youdin, A. Rapid planetesimal formation 
in turbulent circumstellar disks. Nature 448, 1022-1025 
(2007) 

20. Youdin, A. N., Goodman, J., Streaming Instabilities in 
Protoplanetary Disks, Astrophys. J. 620, 459-469 (2005) 

21. Armitage, P. J., Astrophysics of Planet Formation, Cam¬ 
bridge, UK: Cambridge University Press (2010) 

22. Benz, W., Asphaug, E. Catastrophic Disruptions Revis¬ 
ited. Icarus 142, 5-20 (1999) 

23. Leinhardt, Z. M., Stewart, S. T. Full numerical simula¬ 
tions of catastrophic small body collisions. Icarus 199, 
542-559 (2009) 

24. Lissauer, J. J., and 38 colleagues A closely packed sys¬ 
tem of low-mass, low-density planets transiting Kepler-11. 
Nature 470, 53-58 (2011) 

25. Ward, W. R., Protoplanet Migration by Nebula Tides 
Icarus 126, 261-281 (1997) 


10 I www.pnas.org- 


Footline Author 


